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ANALYSIS  OF  ANOMALOUS  RESISTIVITY  DURING  THE  CONDUCTION 
PHASE  OF  THE  PLASMA  EROSION  OPENING  SWITCH 


I.  INTRODUCTION 


The  plasma  erosion  opening  switch  ( PEOS )  is  a  fast  opening 
vacuum  switch  which  has  been  used  to  compress  the  output  pulse 
from  conventional  pulsed  power  generators.1-®  It  is  also  being 
used  in  research  for  developing  high  power  inductive  storage 
pulsed  power  generators. ^  The  PEOS  consists  of  a  low  density 
( n±  -  1013  cm  3)  C++  plasma  injected  (VD  -  10  cm/ps)  between  the 
electrodes  of  a  coaxial  transmission  line  as  shown  schematically 
in  Fig.  1.  The  plasma  is  classically  collisionless  on  the 
timescale  of  interest  (~  100  ns).  However,  the  plasma  has  been 
found  to  conduct  current  (-  1  MA)  across  a  strong  (-  10-50  k  G ) 
se 1 f -gene rated  magnetic  field  until  a  current  threshold  is 
reached  and  then  to  open  (i.e.  cease  to  conduct)  quickly 
(~  10  ns).'  The  current  threshold  is  determined  by  the  injected 
plasma  properties  and  the  geometry.  A  theory  which  describes  the 
mechanisms  involved  in  setting  this  threshold  current  and  which 
describes  the  opening  mechanisms  is  found  elsewhere.^®  Here 
field  penetration  and  current  conduction  in  the  body  of  the 
plasma  during  the  conduction  phase  of  the  PEOS  will  be  addressed. 
This  is  important  for  understanding  PEOS  scaling. 

While  the  PEOS  is  conducting,  electrons  are  emitted  from  the 
cathode  and  enter  the  plasma  after  being  accelerated  across  a 
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cathode  sheath.  Theoretical  analysis  shows  that  the  potential 
drop  across  the  sheath  can  be  large  (-  MV )  .  ^  In  fact,  the 
floating  potential  of  the  plasma  can  be  much  larger  than  the 
induced  voltage  across  the  plasma.  These  energetic  electrons  are 
the  current  carriers  and  are  also  capable  of  driving 
instabilities  in  the  plasma.  Magnetic  field  probe  data1^  shows 
that  the  current  channel  in  a  PEOS  is  many  times  wider  than  the 
collisionless  skin  depth  of  the  plasma  (~  0.1  cm).  in  the 
absence  of  anomalous  resistivity  radial  current  flow  can  only  be 
accomplished  by  Ez  x  B q  drift,  where  Ez  results  from  charge 
separation.  Under  these  conditions  PIC  code  simulations^ , 14 , 
fluid  code  simulations-^-  ^  ^  ®  ^  ^  and  analytic  fluid  work^l  have  all 
predicted  current  channel  widths  that  are  narrower  than  observed 
in  experiments.  PIC  code-^-®  and  fluid  code-^  simulations  with 
anomalous  collisions  modeled  by  introducing  a  collision  term  in 
the  electron  equations  of  motion  have  compared  well  with 
experimental  data  with  an  anomalous  collision  frequency,  va ,  on 
the  order  of  wce- 

Because  the  switch  carries  a  large  current,  the  magnetic 
field  is  large  and  the  electron  gyration  radius  is  small  compared 
with  the  anode-cathode  separation.  The  ions  are  essentially 
unmagneti 2ed .  As  stated  above,  without  anomalous  collisions  due 
to  instabilities  electrons  must  move  with  E  x  B  velocities,  and 
actually  are  tied  to  the  magnetic  field  lines. 

Such  a  motion  can  be  shown  to  lead  to  a  relation  between  ne 
and  B  that  must  be  satisfied  along  electron  trajectories.  This 
constraint  relating  ne  and  B  can  be  combined  with  other  plausible 
assumptions  to  show  that  the  electrons  cannot  cross  the  plasma 
from  the  cathode  to  the  anode  but  must  flow  essentially  on  lines 
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of  constant  r  in  contradiction  to  the  direct  observation  of 
current  flow  across  the  PEOS  plasma.  This  contradiction  then 
provides  the  reductio  ad absurdum  argument  that  motivates  the 
belief  that  there  must  be  a  mechanism  in  the  plasma  to  disturb 
the  pure  E  x  B  motion  of  the  electrons.  In  this  paper  it  is 
suggested  that  the  mechanism  involves  a  collection  of  unstable 
ion  acoustic  waves  that  interacts  with  the  electrons  by  wave 
particle  interactions  to  provide  such  anomalous  collisions.  The 
collisions  produce  the  penetration  of  current  and  magnetic  field 
into  the  plasma  at  an  anomalously  large  rate.  They  explain  the 
ability  of  the  plasma  to  carry  the  large  currents  across  the 
strong  magnetic  field  as  observed. 

The  reductio  ad  absurdum  argument  is  now  presented.  Assume 

that : 

(1)  There  are  no  electron  collisions  and  that  the  electrons 
move  with  the  E  x  B  velocity.  This  is  the  primary  assumption 
whose  absurdity  will  be  shown.  On  the  basis  of  this  assumption 
it  can  be  stated  that  the  magnetic  field  is  frozen  into  the 
electron  fluid.  By  a  standard  proof  it  can  also  be  shown  that 
the  flux  in  a  magnetic  tube  remains  constant  in  time  as  the  tube 
moves  with  the  E  x  B  velocity.  That  is,  for  cylindrical  geometry 

n  r 

— ^ —  =  const.  (la) 

B 

on  the  electron  trajectory  even  for  a  nonsteady  state. 
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(2)  The  plasma  is  charge  neutral.  Usually  the  density  is 
high  enough  for  this  to  be  a  good  assumption.  On  the  basis  of 
this,  Eq.  (la)  becomes 

Zn .  r 

-  =  const.  (lb) 

B 

on  the  electron  trajectory,  where  Z  is  the  ion  charge  state. 

(3)  The  ions  move  very  little  during  the  switch  conduction 
time.  This  can  be  shown  to  be  the  case  for  many  short  conduction 
time  switches  since  the  J  x  B  force  on  the  ions  moves  them  only  a 
short  distance.  Also  their  motion  of  injection  in  the  original 
plasma  is  too  slow  to  move  them  far.  As  a  consequence,  n^  in 

Eq.  (lb)  can  be  taken  to  be  that  of  the  plasma  before  the  power 
pulse  arrives  at  the  switch. 

(4)  The  plasma  is  in  a  quasi-steady  state  so  that  the 
electron  current  lines  are  coincident  with  the  electron 
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trajectories.  From  assumption  (3)  the  ion  current  is  negligible, 
so  that  the  electron  trajectories  are  the  total  current  lines. 

Now  from  cylindrical  geometry  together  with  Ampere's  law  one  can 
conclude  that  B  must  vary  as  1/r  along  current  lines  and  thus 
along  electron  trajectories.  Hence,  Eq.(lb)  becomes 


Zn. r  =  const. 
i 


(  lc  ) 


$ 
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along  electron  trajectories,  where  n^  is  the  initial  ion  density. 

The  argument  can  now  be  completed  as  follows.  If  Zn^  is  a 
function  of  r  alone  and  Zn^  does  not  vary  as  1/r^,  then  the 
electrons  must  move  keeping  r  a  constant.  It  is  more  plausible 
experimentally  that  Zn^  is  independent  of  r  than  varying  as  1/r2, 
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so  that  the  electrons  are  constrained  from  changing  radius  and 
prevented  from  moving  from  cathode  to  anode  in  contradiction  with 
experiment.  Even  if  Zn^  is  not  a  function  of  r  alone,  it  is 
difficult  to  envision  a  reasonable  profile  that  will  allow 
electrons  to  cross  the  plasma  in  light  of  the  constraint 
expressed  in  Eq.  (lc). 

This  argument  requires  certain  comments.  First,  the 
argument  depends  heavily  on  the  cylindrical  geometry. 

Theoretical  analysis  involving  planar  geometry  will  not  encounter 
constraint  (lc).  Second,  assumption  (2)  concerning  charge 
neutrality  is  not  rigorous.  At  lower  densities  it  is  possible  to 
violate  charge  neutrality  over  narrow  channels  (a  few  skin 
depths)  through  which  current  can  pass.  In  this  case  Eq.  (la) 
must  still  be  satisfied  but  Eqs.  (lb)  and  (lc)  need  not  be 
satisfied.  This  appears  to  be  what  is  happening  in  the  PIC  code 
simulations  of  the  PEOS.  However,  for  larger  densities  the 
channel  is  very  narrow  in  contrast  to  the  broad  channels  observed 
in  experiments.  Narrow  current  channels  with  high  current 
density  would  enhance  the  development  of  the  plasma  instabilities 
described  in  this  paper  and  produce  collisions  to  broaden  the 
channel.  These  instabilities  could  be  three  dimensional  (as  is 
the  one  described  herein)  and  would  not  occur  naturally  in  two 
dimensional  simulations. 


The  conclusion  is  that  the  electrons  cannot  carry  out  a  pure 
E  x  B  motion.  There  needs  to  be  an  anomalous  collision  rate,  va , 
that  disturbs  their  motion  so  that  they  are  governed  by 
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For  va  to  make  a  large  enough  difference,  we  see  that  \>a  must  be 
of  order  «ce  -  eB/mec,  the  electron  cyclotron  frequency.  The 
necessary  size  of  va  has  been  confirmed  by  including  va  in 
numerical  calculations  of  E,  B  and  the  fluid  velocity  field,  v, 
which  showed  that  va  must  be  >  .2  wce  to  get  current  patterns 
similar  to  those  observed.  For  typical  PEOS  parameters 
cope  =  ( 4  ne2ne/me )  -  2  x  1011  s_1  and  wce  ~  1-5  x  1011  s-1  at 

the  end  of  the  conduction  phase. 

II.  THE  ION  ACOUSTIC  INSTABILITY  (LINEAR) 

A  possible  mechanism  to  produce  anomalous  collisions  is  ion 
acoustic  waves  driven  unstable  by  the  E  x  B  velocity  of  the 
electrons,  vE  =  cE  x  B/B2,  if  vE  exceeds  several  times 
cs  =  (ZTg/mi)1/2,  the  ion  acoustic  speed.  Here  Te  is  the 
azimuthal  electron  temperature  along  the  magnetic  field  lines. 
Such  waves  are  unstable  for  all  wave  lengths  down  to  the  Debye 
length.  Since  the  wave  lengths  are  in  general  small,  a  local 
theory  can  be  applied  neglecting  any  inhomogeniety  in  the  plasma 
First  the  linearized  theory  of  this  instability  will  be 
cons i de  r ed . 

The  local  E  field  can  be  eliminated  by  transforming  to  a 
frame  moving  with  the  velocity  vE  in  the  E  x  B  direction.  In 
this  frame  the  ions  have  a  mean  drift  velocity  vG  =*  -vE.  The 
ions  are  treated  as  cold  and  the  plasma  is  treated  as 
homogeneous.  A  coordinate  system  is  chosen  with  x  axis  in  the 
direction  of  vE  and  z  axis  in  the  direction  of  B  as  shown  in 
Fig.  2.  (Note  that  the  z  direction  differs  from  the  conventional 


one  along  coaxial  geometry. )  The  electrostatic  wave  dispersion 
relation  for  a  wave  varying  as  exp [ i ( k • r-wt ) ]  is 


e  ( k,  to  )  =  k 
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where  k±  =  kxx  +  kyy,  b  =  k2T±me c2/e2B2 ,  E.  =  c*>' / 1  |  vg  , 

2 

vg  =  Tz/nie  and  Z(x)  is  the  plasma  dispersion  function.  Tx  and 
T ^  (=  Te )  are  the  electron  temperatures  perpendicular  and 
parallel  to  B  respectively.  Taking  b  <<  1  and  E  <<  1,  yields 
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Dropping  the  imaginary  term,  it  is  found  that 
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(  5a  ) 


If  k  XD  <  <  1 ,  then 

<o  =  k*v0  ±  iklcs  (  5b) 

where  AD  *  ve/wpe  is  the  Debye  length  and  cs  =  60p^AD  -  (ZT2/m^ 
is  the  ion  acoustic  speed.  Since  kz  <  0  and  kz  >  0  are 
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equivalent,  k2  >  0  is  chosen  for  simplicity.  Keeping  the 
imaginary  term  in  next  order,  the  growth  rate,  y,  is  found  to  be 
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kcs  £.  exp  (  -1/2) 


(  7a  ) 


where  now 


I  k  •  v  -  k  c 
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In  the  lab  frame  the  unstable  wave  always  propagates  in  that 
direction  along  k  making  an  angle  less  than  90°  with  the  E  x  B 
velocity  vE  (i.e.  -v0).  Since  waves  with  k  and  minus  k  are 
equivalent,  we  discuss  only  those  with  k  •  vD  <  0,  oo '  -  k  •  vQ  =  kcs. 

In  general  j  v0  j  >>  cs  so  that  the  kcs  term  is  negligible 
unless  cost  <<  cs/|v0[,  where  $  is  the  angle  which  k  makes  with 
vE.  Neglecting  the  cs  term,  the  condition  £  -  kxvD/kzve  <  1 
becomes  cost/cos9  <  |vel/|v0|,  where  9  is  the  angle  between  k  and 

B.  The  relation  £,  =  1  determines  a  plane  through  the  vector  E 
and  making  an  angle  tan“^(  |ve[/|v0|)  with  the  E,B  plane.  The 
maximum  growth  rate  ym  =  (n/8)^/^kcs  is  near  this  plane.  If  the 
condition  kXD  <<  1  is  relaxed,  the  maximum  growth  is  found  to  be 
of  order  0.15(*)p^,  corresponding  to  a  typical  growth  time  of 


2.7  ns  for  =  lO15  cm-3  and  C+  .  Although  y  is  considerably 
smaller  than  the  anomalous  collision  frequency  \>a  that  was 
discussed  in  the  introduction,  it  will  be  shown  that  the 
nonlinear  saturation  of  the  wave  is  relatively  weak,  so  that  the 
waves  can  produce  an  anomalous  collision  rate  which  is  faster 
than  their  growth  rate. 


III.  THE  WAVE  PARTICLE  INTERACTION 


The  interaction  of  these  ion  acoustic  waves  will  now  be 
considered  and  it  will  be  demonstrated  that  this  interaction 
leads  to  an  effective  collision  frequency.  In  the  laboratory 
frame 


w  =  oo  +  k-v  =  oo  -  k*v  =  kc 
E  os 


so  that  the  wave  appears  as  an  ordinary  ion-acoustic  wave. 
However,  in  the  E  x  B  frame  it  is  a  negative  energy  wave  driven 
unstable  by  resonant  interaction  with  electrons  moving  in  the 
negative  direction  parallel  to  B.  In  the  laboratory  frame  it 
extracts  energy  by  removing  energy  from  the  E  x  B  motion.  This 
removal  is  the  origin  of  the  anomalous  collision  term  in  Eq .  (2) 

In  addition,  the  electrons  are  heated  in  the  parallel  direction. 
To  see  this,  consider  the  energy,  W^,  and  the  momentum,  P^, , 

of  a  wave  of  amplitude  E  in  the  moving  E  x  B  frame: 
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Thus,  in  the  moving  frame  the  unstable  waves  are  negative 


energy  waves,  propagating  in  the  negative  k  direction,  with 
momentum  parallel  to  k.  In  the  lab  frame 
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(  11  ) 


P  =  P'  ( 12 ) 

w  w 

Thus,  in  the  lab  frame  the  waves  are  positive  energy  waves, 
propagating  in  the  plus  k  direction  with  momentum  in  the  positive 
k  direction. 

Consider  now  the  interchange  of  momentum  and  energy  of  a 
wave  with  the  electrons  and  the  ions  in  E  x  B  frame.  Now  kz  >  0 
and  the  wave  propagation  velocity  along  z  is  k*v0/kz  <  0.  This 
wave  interacts  with  electrons  moving  at  this  same  velocity  in  the 
z  direction  in  the  moving  frame.  The  electrons  are  accelerated 
by  this  wave,  gaining  momentum  and  energy  at  rates 


The  resonant  electrons  gain  a  positive  amount  of  energy  in 
amplifying  the  wave.  They  gain  momentum  in  the  negative  k 
direction.  (These  resonant  electrons  move  in  this  direction.) 

The  ions  are  unaffected  except  through  their  nonresonant 
interaction  with  the  wave. 

In  the  lab  frame  the  electrons  lose  energy  to  the  wave. 

This  arises  by  their  E  x  B  velocity  being  slowed  down,  giving  up 
energy  both  to  the  wave  and  also  to  the  parallel  part  of  the 
electron  energy.  There  is  consequently  a  force  on  the  electrons. 
The  part  of  this  force  perpendicular  to  B  is 


Fx  =  -12  n 


.  2 

k  ,  a) 
x  Pi 


.22  k  v 
k  c  s  z  e 


In  general,  the  waves  are  symmetrical  with  respect  to  the  vE,B 
plane  so  that  any  component  of  Fx  perpendicular  to  v£  cancels 
out.  The  component  of  Fx  along  v£  is  in  the  opposite  direction 
to  the  E  x  B  motion.  Denoting  this  component  by  -vanemevE,  and 
setting  E  =  -ik<J>  yields 


.2  ~  o 
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The  heating  in  the  parallel  direction,  Eq .  (13)  can  be  reduced  to 


the  form 


dw '  2 

■■  ,-e  =  -F,  •  v_  =  v  n  m  v_ 
dt  -1  E  a  e  e  E 


where  Eq.  (15)  has  been  employed.  This  corresponds  to  the 
frictional  heating  of  the  electrons  by  the  anomalous  collisions 
with  the  ions. 

The  wave-particle  interaction  in  the  laboratory  frame  is 

illustrated  in  Fig.  3.  The  wavevector  k  is  chosen  so  that 

k z  >  0.  The  ion  acoustic  wave  propagates  along  k  with  velocity 

cs.  The  condition  for  an  electron  with  guiding  center  velocity 

u  =  v_  +  v  z  to  be  in  resonance  with  the  wave  is 
E  z 

k,*v_+kv=w=kc.  (18) 

-1-  E  2  z  s 

Neglecting  cg  compared  to  v E ,  we  have 


so  that  the  resonant  electron  moves  in  the  negative  z  direction. 
The  vector  u  makes  an  angle  of  nearly  90°  with  k  (if  k  is  in  the 
vE , B  plane ) . 

Because  of  the  slope  of  the  electron  distribution  function 
in  the  parallel  direction,  the  resonant  electrons  extract 
negative  momentum  from  the  wave  (see  Fig.  3).  Since  the  momentum 
of  the  wave  is  in  the  positive  2  direction  this  means  that  the 
wave  must  amplify,  the  energy  coming  from  the  vE  motion  cf  the 
electrons.  (Figure  3  has  been  drawn  for  k  in  the  vE,B  plane  but 
k  can  be  in  a  general  direction  and  amplification  will  occur 
provided  that  k*vE  >  0.) 

From  Eq .  (16)  it  can  be  seen  that  if  kx  ~  k2  -  1/XD  and 
e$/Tz  =  1,  then  \>a  =  wpe '  a  value  typically  of  order  of  the 


electron  cyclotron  frequency  as  required  by  the  discussion  in  the 
introductory  section.  The  magnitude  of  e$/Tz  is  estimated  in 
the  next  section. 


IV.  NONLINEAR  THEORY 


In  this  section  the  nonlinear  amplitudes  of  the  unstable  ion 
acoustic  waves  are  estimated  in  order  to  determine  the  actual 
magnitude  of  va  from  Eq .  (16).  Because  of  the  negative  energy 

character  of  the  waves  they  do  not  saturate  by  quasilinear 
saturation.  Instead  the  resonant  particles  are  accelerated  and 
heated  in  the  direction  parallel  to  B  as  the  waves  are  amplified. 
When  the  waves  reach  a  certain  amplitude,  they  are  saturated  by 
induced  scattering  off  of  ions.  The  equation  for  the  evolution 
of  the  wave  amplitudes  is  given  by  Kadomtsev^0  for  the  limit 
Ti  <<  Te  =  Tz  and  can  be  written  as 
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where  I^d^k  is  the  mean  contribution  to  <$*->  from  waves  with  k  in 
d3k,  k"  -  k  -  k',  w"  *  is  the  frequency  cf  the  wave 

with  wave  vector  k  and  6'  is  the  derivative  of  the  delta 
function.  The  right  hand  side  represents  the  rate  of  transfer  of 
intensity  from  the  k'  mode  to  the  k  mode  by  induced  scattering. 

In  a  steady  state,  when  modes  are  saturated,  energy  should  be 
scattered  predominantly  from  unstable  modes  to  damped  modes, 
although  there  should  be  some  scatter  among  the  unstable  mode e  as 
well  as  among  the  damped  modes.  The  unstatle  modes  are  sera:  ate;: 


from  the  damped  modes  by  the  plane  in  k  space  through  the  origin 
and  perpendicular  to  vE. 

Making  use  of  the  result  that  =  kcs  and  from  the  delta 
function  factor  which  shows  that  k  =  k',  we  may  reduce  Eq .  (20) 
to  the  form 
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where  p"  =  k*k'/kk'  and  where  we  have  introduced 


(21  ) 


(22) 


In  terms  of  which  is  dimensionless,  the  expression  for  va  may 

be  written  as 


v 
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k2« 

d3k 

lkJk 


(23) 


Equation  (21)  is  an  i n te g r o-d i f f e r en t i a  1  equation  for  4>k , 
which  must  be  solved  numerically.  However,  its  interpretation 
is  clear.  Waves  with  positive  kx  are  unstable  with  yj<  >  0,  so 
that  their  amplitude  grows  due  to  the  y^  term  in  Eq .  (21). 

Similiarly,  y^  <  0  for  the  damped  modes  with  kx  <  0,  and  their 
amplitude  decreases  due  to  the  y^  term.  The  nonlinear  scattering 
term  on  the  righthand  side  of  Eq .  (21)  transfers  energy  between 
all  the  modes  with  k'  =  k.  In  equilibrium  the  transfer  of  energy 


from  unstable  to  stable  modes  must  balance  the  growth  of  the 


unstable  modes,  and  the  damping  of  the  stable  modes.  The  sign  of 
the  right  hand  side  must  be  such  as  to  give  this  balance.  The 
transfer  of  energy  among  the  unstable  modes  can  be  regarded  as  a 
rearrangement  of  the  energy  of  the  unstable  modes,  but  does  not 
contribute  to  the  equilibrium  balance. 

To  get  an  idea  as  to  the  magnitude  of  the  energy  in  the 
unstable  modes,  consider  the  assumption  that  ^  is  independent  of 
angle  for  the  unstable  modes  (kx  <  0)  and  a  function  of  k  alone, 
say  Similiarly,  assume  that  4^  is  a  different  function  of 

k  alone,  say  4>s(k),  for  the  stable  modes.  This  assumption  for 
will  not  solve  the  full  Eq .  (21)  of  course,  but  it  can  be  chosen 
to  make  two  moments  of  this  equation  vanish.  We  average  Eq .  (21) 
for  this  choice  for  over  all  unstable  directions  kx  >  0,  and 

again  over  all  stable  directions  kx  <  0,  to  obtain  two 
differential  equations  for  4>u  and  4>s .  [This  corresponds  to  the 
Eddington  approximation  in  the  theory  of  radiative  transfer.  It 
can  be  generalized  by  expanding  ♦  in  a  series  of  Legendre 
functions  and  taking  Legendre  moments  of  Eq .  (21).] 

These  assumptions  can  be  written  as 

*.  =  *U(k)  k  >  0  ,  (24a) 

k  x 

$.=«>*(  k  )  k  <  0  (  24b) 

k  x 

In  caking  the  angular  averages  we  must  neglect  the  coupling  of 
two  unstable  modes,  or  two  damped  modes,  since  physically  this 
coupling  does  not  lead  to  a  change  in  either  the  stable  or 
unstable  energy.  This  fact  is  not  represented  by  our 
approximation. 


Substituting  Eqs.  (24a)  and  (24b)  into  Eq .  (21)  and 


integrating  first  over  the  hemisphere  kx  >  0  with  the  angular 
integral  over  k'  restricted  to  kx'  <  0  yields 
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Carrying  out  the  two  integrals  over  the  opposite  hemispheres 
yields 
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Here  <Yk>u  and  <rjc>s  are  the  growth  rates  given  by  Eq .  (7) 

angularly  averaged  over  the  unstable  and  stable  hemispheres. 
Their  averages  are  given  by 


<Vu 


,  ^  <£exp(-£,2/2)>^kcs 
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=  0.184  f ( a)  kcs  , 
<Yjc>s  =  -  0.184  f  (  a )  kcs 


(  26b) 


where  a  *  v^/ve,  E^(x)  is  the  exponential  integral^  and  f(a)  is 
normalized  to  unity  at  a  -  1.  For  small  a,  f(a)  -  2 . 1 7 aln ( 1/a ) , 
while  for  large  a,  f(o)  -  2.17/a. 

Assuming  a  steady  state  in  Eqs.  (25)  yields 
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where  0  =  0 . 22 ( ZTe/T^ ) f ( vE/ve  )  .  The  solutions  of  these  equations 
are 

*S  =  3  In ( k/k  .  )  ,  ( 28a  ) 


$U  =  6  ln(kD/k) 


(  28b  ) 


Boundary  conditions  have  been  chosen  so  that  *s  is  zero  at 
some  minimum  value  of  k  denoted  by  km^n.  That  km^n  is  not  zero 
is  a  reflection  of  the  fact  that  for  small  k  the  growth  rate  is 
too  slow  to  build  up  the  waves  instantaneously  to  the  steady 
state  value.  (km^n  should  decrease  in  time  as  1/t . )  Also  the 
condition  that  4>u  *  0  at  k<=  kD  =  1/XD  has  been  employed  since  at 
shorter  wavelengths,  the  waves  are  damped  by  ion  Landau  damping. 
However,  a  finite  amplitude  for  the  unstable  waves  can  be 
sustained  by  nonlinear  coupling. 

One  may  now  substitute  the  expressions  for  Eqs.  (24), 

into  Eq .  (23)  to  evaluate  va  .  There  results 
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where  <kx^/|kz|k>  represents  the  angular  average  of  kx^/()k2|k). 
This  average  is  logarithmically  divergent  for  small  kz.  However, 
because  the  growth  rate,  given  in  Eq.  (7),  is  exponentially  small 


for  small  k 


z,  we  may  reasonably  cut  this  average  off,  say  at 
kz  =  0.1.  This  yields  <kx^/|kz|k>  =  0.9.  Substituting 
Eqs.  (28a)  and  (28b)  into  Eq .  (29),  we  find 
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(  30  ) 


As  noted  above  km^n  decreases  with  time,  so  va  should  slowly 
increase  with  time. 

Although  the  procedure  employed  to  solve  Eq.  (21)  is 
admittedly  very  approximate,  it  does  bring  out  the  physics.  The 
resulting  value  for  the  anomalous  collision  frequency,  va ,  is 
correspondingly  approximate,  but  it  should  be  correct  in  order  of 
magnitude.  Equation  (30)  indicates  that  the  ion  acoustic  mode 
driven  unstable  by  the  large  E  x  B  drift  can  reach  an  amplitude 
that  is  sufficiently  large  to  produce  a  va  as  large  as  the 
electron  plasma  frequency,  «pe ,  even  though  the  growth  rate  of 
the  mode  is  of  order  cjp^  .  The  reason  that  this  is  possible  is 
that  the  mode  coupling  that  saturates  the  waves  is  so  weak  that 
the  waves  grow  to  a  very  large  amplitude  independent  of  its 
linear  growth  rate.  Further,  because  the  waves  propagate  at  the 
ion  acoustic  speed,  they  can  achieve  this  large  amplitude  without 
leaving  the  plasma. 

It  should  be  noted  that  it  was  assumed  that  the  mode  is 
strongly  unstable,  so  that  vE  >>  cs  and  so  that  the  c5  term  oar. 
be  neglected  in  Eq.  (7b).  If  \>a  is  too  large,  v E  will  be 
reduced,  and  cs  will  increase  through  heating  of  the  electron 
temperature  parallel  to  B.  Eventually  the  unstable  angular 
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region  in  the  kx  >  0  hemisphere  will  be  reduced  and  va  will  be 
correspondingly  reduced. 

Finally,  it  should  be  noted  that  the  equations  employed  in 
the  derivation  of  the  saturated  mode  amplitude  and  of  va  are 
based  on  the  small  amplitude  limit  e$/Te  <<  1.  If  e$/Te  ~  1, 
many  more  nonlinear  effects  must  be  included.  Hence,  the 
coefficient  of  o>pe  is  not  accurate.  At  most  it  can  be  said  that 
e?  -  Te  and  that  \>a  -  «pe . 

V.  CONCLUSIONS 

It  has  been  shown  that  an  anomalously  large  electron-ion 
collision  frequency  is  needed  to  allow  the  flow  of  current  across 
the  large  magnetic  field  in  the  PEOS .  From  general 
considerations  this  collision  frequency  is  shown  to  be  of  order 
of  a  few  tenths  of  the  electron  cyclotron  frequency. 

It  is  proposed  here  that  one  mechanism  for  producing  this 
anomalous  collision  frequency  is  a  large  amplitude  ion  acoustic 
waves.  These  waves,  when  propagating  in  a  direction  which  makes 
an  angle  of  less  than  90  degrees  with  respect  to  the  E  x  B 
direction,  are  driven  unstable  by  the  large  E  x  B  velocity  of  the 
electrons  through  the  ions.  The  growth  rates  of  these  waves  is 
comparable  with  the  ion  plasma  frequency.  In  most  cases  this 
frequency  is  very  large  compared  with  the  time  scales  for  switch 
operation,  so  that  there  is  time  for  many  exponentiations  of  the 
wave  amplitude.  Since  the  waves  propagate  at  the  ion  acoustic 
.-peed,  they  do  not  propagate  very  far  during  their  growth  times 
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so  that  they  are  absolutely  rather  than  convectively  unstable. 

Examination  of  the  quasilinear  wave-particle  interactions  of 
these  waves  has  shown  that  they  remove  momentum  and  energy  from 
the  directed  motion  of  the  electrons  perpendicular  to  the 
magnetic  field  at  a  rate  that  can  be  represented  by  the  anomalous 
collision  frequency  va  .  The  resulting  frictional  energy  goes  to 
heat  the  thermal  motion  of  the  electrons  parallel  to  the  magnetic 
field.  The  thermal  motion  perpendicular  to  the  field  is 
uncnonged  because  of  the  adiabatic  invariance  of  the  magnetic 
moment  of  the  electrons.  When  the  modes  damp,  they  transfer  the 
momentum  to  the  ions. 

To  determine  the  magnitude  of  va  it  is  necessary  to 
determine  the  saturated  amplitude  of  the  waves.  There  is  no 
saturation  of  the  waves  by  quasilinear  wave-particle  interactions 
because  the  energy  of  the  waves  is  effectively  negative.  The 
only  quasilinear  saturation  effect  is  to  heat  the  parallel  motion 
of  the  electrons  so  as  to  raise  the  acoustic  speed  cs.  This 
speed  is  usually  much  smaller  than  the  driving  E  x  B  velocity  v  j- . 
The  saturation  mechanism  is  consequently  the  much  weaker  one  of 
nonlinear  scattering  of  the  unstable  ion  waves  into  stably 
propagating  directions.  The  nonlinear  i nteg r od i f f e r en t i a  1 
equation  which  determines  the  saturation  was  approximated  by  a 
simple  Eddington  type  assumption,  in  order  to  derive  an 
approximate  value  for  the  saturated  amplitudes. 

The  result  is  that  because  the  saturation  process  is  weak, 
the  ion  acoustic  waves  grow  to  large  amplitudes  (e$/Te  -  1',  and 
the  resulting  anomalous  collision  frequency  is  of  order  of  the 
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electron  plasma  frequency.  In  most  cases  of  PEOS  operation  this 
is  of  the  same  order  of  magnitude  as  the  electron  cyclotron 
frequency  and  is  large  enough  to  explain  the  cross  field  current 
It  is  thus  concluded  that  the  unstable  ion  acoustic  waves 
are  strong  enough  to  produce  the  collisions  necessary  to  explain 
the  cross  field  currents  and  the  magnetic  field  penetration  into 
the  PEOS  plasma. 
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STORAGE 

INDUCTOR 


Fiq.  i  -  The  geometry  of  the  PEOS .  The  directions  of  the 
electric  field  and  the  resulting  E  x  B  drift  at  a  point  are 
indicated.  The  coordinate  system  employed  in  Secs.  I I-V  is 
d : si  laved. 
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Fig.  2  -  The  orientation  of  the  coordinate  system  relative  to  E, 
B  and  vE  is  shown.  Also  for  a  general  wavevector  k,  the  various 
components  are  given.  k  is  chosen  in  the  unstable  direction. 
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